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Within a variational approach to solve the Gross-Pitaevskii equation we investigate dynamical 
properties of a rotating Bose-Einstein condensate which is confined in an anharmonic trap. In 
particular, we calculate the eigenfrequencies of low-energy excitations out of the equilibrium state 
and the aspect ratio of the condensate widths during the free expansion. 

PACS numbers: 03.75.Kk 



I. INTRODUCTION 

The rotation of a Bose-Einstein condensate (BEC) has 
raised enormous interest, both experimentally and theo- 
retically over the last few years [lj . The BEC is a super- 
fluid, hence its rotational behavior is quite different from 
that of a normal fluid. Its properties strongly depend on 
the specific trap in which it is confined. In harmonic ro- 
tating traps centrifugal forces may break the confinement 
when the rotation frequency gets close to the trapping 
frequency. To overcome this centrifugal barrier, one has 
to increase the strength of the radial confinement. The 
most natural approach relies on adding a radial quar- 
tic potential to the harmonic trap. Such an anharmonic 
trap was realized in a recent experiment in Paris at the 
Ecole Normale Superieure (ENS) in the group of J. Dal- 
ibard, where the centrifugal limit can be exceeded slightly 
0, This combined trap changes both the thermo- 
dynamical and the dynamical properties of the conden- 
sate drastically, in particular in the fast rotation regime. 
Nevertheless, the analysis of the experimental data has 
been mainly performed by considering only a harmonic 
trap. This approximation is valid for a slow rotation fre- 
quency, but a correct interpretation of the experimental 
data in the fast rotation regime necessitates an analysis 
for the full anharmonic trap. Whereas the proceeding pa- 
per [J] elaborates within the semiclassical approximation 
the thermodynamical properties of such a rotating BEC 
as, for instance, the critical temperature and the heat 
capacity, the present work investigates the correspond- 
ing dynamical properties. To this end we proceed as 
follows. Section [II] introduces the basic Gross-Pitaevskii 
equation for the condensate wave function of a rotating 
trapped BEC. Neglecting the emergence of vortices, we 
determine in Section IIIII the stationary equilibrium den- 
sity of a rotating condensate within the Thomas-Fermi 
approximation. In Section IIVI we then use a variational 
approach to investigate the hydrodynamic modes of the 
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condensate and its free expansion after switching off the 
trap. All results are obtained for the specific parameters 
of the Paris experiment @, Q • 



II. HYDRODYNAMIC FORMALISM 

At zero temperature, the dynamics of a BEC is de- 
scribed by the Gross-Pitaevskii Lagrangian density [f| 

dt 2M 



-F(x)V>*V> ' 



2irTi 2 a s 
M 



(1) 



Here, V{x) denotes the trapping potential, a s the s-wave 
scattering length, and M the atomic mass. The transfor- 
mation into the rotating frame has the consequence that 
an angular momentum — 1~2 L is added to the Lagrangian 
([I]) , where the rotation vector is denoted by = f2e z and 
L represents the one-particle quantum mechanical angu- 
lar momentum operator 0, [|[ . Thus, the Euler-Lagrange 
equation yields the Gross-Pitaevskii (GP) equation 
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Performing the decomposition tp = \ip\e lS we work out 
the hydrodynamic aspects of the system, with condensate 
density n = \ip\ 2 and superfluid velocity v s = h'VS/M. 
Basic manipulations yield for them the coupled hydrody- 
namic equations which consist of a continuity equation 
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and the Eulcr equation 
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FIG. 1: Chemical potential versus rotation fre- 
quency for the values of the Paris experiment @, 0| . 
The solid line corresponds to (|10[) . the short-dashed 
and the long-dashed line to the approximation (|13p 
and (H5), respectively. 



FIG. 2: Aspect ratio of condensate radii versus ro- 
tation frequency for the values of the Paris experi- 
ment 2, 3]. Solid line corresponds to the Thomas- 
Fermi radii R^ F /R^ , dashed line to the stationary 
radi i Wq z f /WqJ determined from Eq. J2TJ in Sec- 
tion Hvll 



where 7(x,0) = V{x) - Mfl 2 (x 2 + y 2 )/2 denotes the 
effective potential in the rotating frame. Here, the solid- 
body velocity is denoted by v s b = xx In the following, 
we consider the cylindrical anharmonic trap used in the 
Paris experiment @, [|[ 
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where M is the atomic mass of 87 Rb, u>± — 2ir ■ 
64.8 Hz, uj z = 2 tt ■ 11.0 Hz, k = 2.6 ■ 1Q" 11 Jm" 4 , and 
r± = \J x 2 + y 2 . An amount of 3 • 10 5 atoms has been 
condensed and the highest realized rotation speed was 
O max » 1.04- ux. 

The superfluid is irrotational V x v s = 0, except at sin- 
gularities, whereas V x v s b = — 2f2. Thus the difference 
of both v s — v s b cannot vanish identically. However, the 
superfluid can mimic solid-body rotation by forming an 
uniform array of rectilinear vortices parallel to the axis of 
rotation 0] , which minimizes the difference v s — v s b . In 
the detailed analysis of Ref. 10] it has been shown that 
the triangular lattice is the energetically most favorable 
one. For a dense vortex lattice it is even justified to ne- 
glect in the Euler equation ^ the term proportional to 
the difference v s — v s b 11]. In the following we assume 
that this approximation is valid and focus on investigat- 
ing the effect of the centrifugal term in the potential ([5]). 



III. STATIC PROPERTIES 



Furthermore, we perform the Thomas-Fermi (TF) ap- 
proximation [l3|, i.e. the kinetic energy is neglected in 
©, so that the condensate density reads 
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where the Heaviside function ensures its positivity. 
The chemical potential is determined by the normaliza- 
tion condition N = J d 3 a;n TF . The Thomas-Fermi radii 
for the cylinder symmetric trap potential (O read 
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We work out the spatial integration in cylinder coordi- 
nates, yielding 



N 



2V2M/X 2 f M(u 2 ± -Q 2 ) 

2 J 



where the function 

f( x ) = T^( 1 + x2 ) 2 l--arcsin^= 
16 [ n VI 

5x 3x 3 



(10) 
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We start with analyzing the static properties of such 
a condensate, where we follow Ref. [12} • To this end 
we separate the time-dependence according to n = n(x) 
and S = —/J,t/h+ s(x), where /j, denotes the chemical 
potential. With this we obtain from the Euler equation 
the stationary hydrodynamic equation 
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has the asymptotic properties 

' x/5 - x 3 /35 for x > 1 

3tt(x 2 + 1) 2 /16- x - x 3 for |x| < 1 
I 3tt(x 2 + l) 2 /8 - x/5 for x < -1 . 

When the rotation frequency f2 varies between and 
^max = 1.04- the argument x — M{lo\ — il 2 ) / (2^/kJl) 
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of the function (|TT|) ranges between — oo < x < 3.3. A 
slow rotation corresponds to x 3> 1, where the chemical 
potential reads 
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For the fast rotation regime ~ 
second expansion (fl2|) . where |x| 
the chemical potential 



we have to use the 
<C 1, so we obtain for 
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For an arbitrary rotation frequency ft the chemical poten- 
tial follows directly from numerically inverting (jlOp . Ac- 
cording to Figure [1] we observe that the chemical poten- 
tial vanishes for the experimental data of the Paris exper- 
iment at a critical rotation frequency fl c » 1.07 -uj±. 
Our assumption, that this indicates the emergence of an 
instability, is confirmed by evaluating the TF radii l[5]). 
© and the TF density Q. Figure H] depicts that the TF 
radius in z-direction vanishes with the chemical potential 
as it is proportional to the square root of it, while the 
other radius in the perpendicular plane remains finite. 
Note that the vanishing of one TF radius also signals 
the breakdown of the validity of the TF approximation. 
Thus, our TF value f2 c » 1.07-u/j_ for the critical rotation 
frequency should be modified by a full quantum mechan- 
ical calculation. In Figure [3] the TF density is shown for 
three particular rotation frequencies. For an overcriti- 
cal rotation f2 > u>±, we notice a pronounced dell in the 
density profile. Thus, we can conclude from Figures [TH3] 
that our TF analysis for the stationary condensate den- 
sity qualitatively agrees with the experimental finding 
that the rotation frequency can only be increased until 
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FIG. 3: Thomas-Fermi density (0 in the xw-plane 
for the values of the Paris experiment 0, Hf an d 
various rotation frequencies Q. From top to bot- 
tom the rotation frequencies are Sl/cu± — (solid), 
Q/lu± = 1 (short-dashed) and Q/lj± = 1.06 (long- 
dashed) . 



first order of 6n and Sv we obtain the eigenmode equa- 
tion by taking the time derivative of the Euler equation 
(|4"|) and inserting the continuity equation §3§ 

Mto 2 Sn = -\fi- V(x, fi)] V 2 Sn + W(x, fi.) • VSn , 

, , . (15) 

where the chemical potential is a constant which is de- 
termined from normalizing the stationary density as de- 
scribed in the previous section. The eigenmode equa- 
tion (Ti~5|) has been solved analytically for a harmonic 
trap which is isotropic [l4| or rotationally symmetric [l5| . 
However, investigating Eq. (|15[) for the anharmonic trap, 
we observe that the differential operator does not con- 
serve powers of the radial coordinate r± , so that we may 
not expect simple analytical solutions. In the special case 
of an overcritical rotation fi > ui± an analytical solution 
of (14) was found in Ref. (l7j . 

Another way to obtain the eigenmodes at least approx- 
imatively is a variational approach which is based on 
Ritz's method. This method was already successfully ap- 
plied for a harmonic trap [T(| and is in agreement with 
the solution of (|15p for a large two-particle (S-interaction, 
i.e. the regime where the TF approximation is valid. For 
a simplified one- and two-dimensional case this method 
was already applied for the anharmonic trap ((SJ) [lH [l!^ ] . 
However, a correct description involves a full three di- 
mensional treatment which we present now. 



IV. DYNAMIC PROPERTIES 



Variational Method 



In order to investigate the dynamical properties of the 
condensate, we return to the Euler equation ^ in the 
TF approximation and consider collective excitations of 
the stationary state. Thus, we look for solutions of the 
type n(x, t) — n TF (x) + Sn(x) e luJt , where n TF (x) is the 
stationary TF density ([7]) and 5n(x) e lut a small time- 
dependent deviation. Analogously, the superfluid veloc- 
ity field is decomposed according to v s = v s b + Sv. In 



For our variational approach we follow Ref. [16] and 
use a Gaussian shape of the condensate 
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FIG. 4: Three eigenmodes of the non-rotating con- 
densate (fi = 0) confined in the anharmonic trap 
(JSJ): a) breathing mode, b) and c) scissor modes. 
The solid lines correspond to the solution of the 
variational procedure (I20p . the long-dashed lines to 
its TF approximation (|22[) . and the short-dashed 
lines to the harmonic trap (k = 0). The conden- 
sate clouds visualize the temporal behaviour of the 
respective eigenmodes. 



where W XlVlZ and S XjVtZ denote time-dependent varia- 
tional parameters for the widths and the phases, respec- 
tively. Note that (JTSJ) can only be used in the undercrit- 
ical regime < luj_. For il > lu±_ a central hole emerges 
in the condensate (see Figure [3]) , which is not described 
by the Gaussian trial function. For the ansatz (fi"6|) we 
determine the Lagrangian L = J d 3 x C with the den- 
sity ([T]) and the potential ([5]). Thus, the resulting action 
A = J dt L represents a functional of the trial functions 
W x , y , z and S x ,y. z which are determined according to the 
Hamilton principle. As we find for the phases the rela- 
tion S XjVtZ = —W x , y ^ z /2W XtViZ , we can eliminate them 
and find for the widths the coupled equations of motion 
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where we have introduced the dimensionless variables 
A = uj±/uj z » 6, i] = 1 — (n/oj±) 2 , and k = kfiuj z ja\ 
0.4. The length s cale is set by the harmonic oscilla- 
tor length a z = ^Ti/Muj z , so that W x>y , z -> W X}VtZ a z , 
and the time is scaled by t — > t/u z . The resulting di- 
mensionless interaction parameter is denoted by P = 
yj2/i:Na s /a z . The term proportional to k is due to the 
anharmonicity of the trap and represents the main differ- 
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FIG. 5: Eigenfrequencies of the condensate con- 
fined in the anharmonic trap ([5]) in the TF approx- 
imation dm. 



ence in comparison with the treatment of the harmonic 
trap in Ref. [16| . The Euler-Lagrange equations (fT7)l are 
of the form W x , v , z = -dV e s(W x ,W y ,W z )/dW x , VlX and 
can therefore be regarded as the motion of a fictitious 
point particle in the effective potential 
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where the vector of the width is denoted by W 

(W X ,Wy,W Z ). 



B. Low-Energy Excitations 

At first, we determine the optimal widths of the equi- 
librium state around which the density fluctuates. The 
stationary point is obtained by setting Wrjfc = 0, yielding 
a set of algebraic equations. Because of the rotational 
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tionary widths follow from 
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Considering small deviations of the equilibrium state, we 
perform a Taylor approximation of the effective potential 
V eS (W) sa V cS (W )+5WH5W T /2, where W denotes 
the stationary point and H is the Hessian matrix. The 
square-root of the eigenvalues of H are the low excitation 
frequencies, for which we find the analytical expressions 
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FIG. 6: Aspect ratio of the non-rotating conden- 
sate (fi = 0) for the values of the Paris experiment 
0, 0|- The solid line corresponds to k = 0.4, the 
long-dashed lines to k = 0, and the short-dashed 
lines to the interaction free Bose gas (a s — 0) with 
K = 0.4. 
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where 



1/2 



[X 2 tj + 6kW$ x - 1 



-^2,3] + 8P2.3-P4 



we have introduced the notation P^j 
P / (4:Wq x W^ z ) . The corresponding TF approximation of 
these modes follows from neglecting the l/W 3 -teims in 
(I17p which correspond to the kinetic energy and is here 
related to a large interaction parameter P — > 00. Thus, 
the algebraic equations (flTj|) for the stationary widths re- 
duce to 
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The resulting ratio Wj z F /W^J agrees qualitatively with 
the aspect ratio of TF condensate radii (see Figure O . 
The oscillation frequencies (f2T))) simplify in the TF ap- 
proximation to 



TF 
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Figure |4] shows that the three eigenfrequencies ([20]) and 
their TF approximation Ij22|) for a non-rotating conden- 
sate (fi = 0) increase with the s-wave scattering length 
a s . This effect is due to the presence of the trap an- 
harmonicity, as the eigenfrequencies in a harmonic trap 
(k = 0) tend to a constant in the limit a s — > 00. For 87 Rb 
the s-wave scattering length is about a s = 5 nm, so we 
can conclude from Figure U that the TF approximation is 
sufficient for the Paris experiment [1, [1]. Therefore, the 
resulting strong dependence of the eigenmodes on the ro- 
tation frequency f2 in Figure [5] is only depicted for the 
TF approximated frequencies 




FIG. 7: Aspect ratio of the rotating condensate 
confined in the Paris trap 0, 0| for various ro- 
tation frequencies. The rotation parameter r\ — 
1 — (Q/ojjJ 2 is varied in steps of width 0.2 from 
(lowest line) to 1 (highest line). 



C. Free Expansion 

The equations of motion (fT7|) determine also the ex- 
pansion dynamics of the condensate after switching off 
the trap. Following Ref. [2fjj . we rewrite them within the 
TF approximation, yielding 



1 



X 2 j 1 + AkW^ 2 



b\b z 



(23) 
(24) 



where we have introduced the scaling parameters b± (t) = 
Wr(t)/Wl F and b z (t) = Wj F {t)/Wj F . The ini- 
tial conditions for the scaling parameters are given by 
6 ± (0) = 6,(0) = 1 and b±(Q) = b z (0) = 0. The ini- 
tial velocities vanish at t = because of the equilibrium 
condition. Since we have A 2 7; + AkW^J 2 > 1 for all 
experimentally realizable rotation frequencies f2, the ex- 
pansion in z-direction is much slower than in radial direc- 
tion. Therefore, we may decompose the scaling parame- 
ter b z = 1 + e and consider e as a small time-dependent 
quantity. In zeroth order in e we obtain from (|23[) for the 
radial scaling parameter 
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With this solution we determine the scaling parameter 
for the expansion in z-direction by integrating e = l/b± 
either numerically or analytically in the limit of slow and 
fast rotation frequencies. With this we approximately 
obtain for Q -C lu z 

t arctan f y/X 2 r] + 4 K Wg k ' 2 t) 
b z « 1 + , — '- 

\n[l + {X 2 r 1 + AKW^ 2 )t 2 } Vtt/uo z 
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and, correspondingly, for f2 « uj z 

b z « 1 + In ^cosh . (27) 

Note the time of flight of the Paris experiment @, [|[ lasts 
18 ms. Figure [6] shows that the anharmonicity k has 
only a marginal effect on the aspect ratio of the widths 
Wj F /Wj F = ^Wg r J F /6 z Wg I 2 F for a non-rotating conden- 
sate (S7 = 0). However, this changes drastically once the 
condensate is rotating. From Figure [7] we read off that 
the effect of the trap anharmonicity upon the aspect ra- 
tio is significantly enhanced in the fast rotating regime. 
Therefore, this effect must be taken into account for any 
analysis of free expansion data for large rotation speeds. 

V. CONCLUSIONS 

In this paper we have analyzed dynamical properties 
of a rotating BEC for the specific parameters of the 
Paris experiment 0, Hj]. With this we have provided 
useful theoretical results for future BEC experiments in 
such anharmonic traps. 

At first, we have determined how the stationary TF 
density profile of the condensate varies with increasing 
rotation frequency Q. In the overcritical rotation regime 
we have found that the density in the center is reduced. 



Thus, the rotation frequency has a theoretical upper 
bound of Q c « 1.07 • lu±, which qualitatively agrees with 
the highest experimentally realized rotation speed of 
!] max w 1.04 • u>± 0, ||. This quantitative agreement is 
astonishing insofar, as we have neglected the presence of 
vortices in our TF analysis. 

Furthermore, we have used a variational approach to 
find the frequencies of the low-energy excitations which 
strongly depend on the rotation frequency Q. The an- 
harmonicity influences the eigenmodes in a characteristic 
way as it enters the equilibrium state of the condensate. 
Finally, we have determined the time dependence of 
the aspect ratio of the condensate widths for the ex- 
pansion out of the trap which is an important quantity 
for analyzing the experimental data. The rotation 
increases significantly the velocity of the expansion in 
the a;y-plane. This effect could be used to determine the 
rotation speed in the experiment. 
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